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Abstract 



We compute the one-loop Casimir energy of gravity and matter fields, obeying various 
boundary conditions, in 5-dimensional S 1 /Z 2 and 6-dimensional T 2 jZy. orbifolds. We discuss 
the role of the Casimir energy in possible radius stabilization mechanisms and show that the 
presence of massive as well as massless fields can lead to minima with zero cosmological constant. 
In the 5-d orbifold, we also consider the case where kinetic terms localized at the fixed points 
are not small. We take into account their contribution to the Casimir energy and show that 
localized kinetic terms can also provide a mechanism for radius stabilization. We apply our 
results to a recently proposed 5-dimensional supersymmetric model of electroweak symmetry 
breaking and show that the Casimir energy with the minimal matter content is repulsive. 
Stabilizing the radius with zero cosmological constant requires, in this context, adding positive 
bulk cosmological constant and negative brane-tension counterterms. 



1 Address after July 1, 2001: Department of Physics, University of Toronto, 60 St George St., Toronto, ON 
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1 Introduction and summary. 



The idea that there are extra spatial dimensions in which gravity and, possibly, some or all 
matter fields can propagate has been the subject of renewed interest in the last few years. The 
Kaluza-Klein and "braneworld" ideas combined have been used to reformulate and address 
every conceivable problem of elementary particle physics |1J - [12]. 

One of the main theoretical issues in theories with extra dimensions is that of determining 
their size. In the absence of a stabilization mechanism, the Casimir energy tends to either 
inflate or contract the extra dimensions, as has been known since the work of refs. [13|, [14| . The 
Casimir energy of various braneworld compactifications has also received attention recently 



T|, [13, |T|, |T|, |0|. 



In this paper, we focus on the Casimir energy and its role in radius stabilization in a partic- 
ular class of Kaluza-Klein compactifications not discussed until recently — field-theory orbifold 
compactifications.]^] Orbifolds are very useful tools for projecting out unwanted massless modes 
and/or breaking (super) symmetries. Several interesting phenomenological models using field 
theory orbifolds were proposed recently |22, 23 1 . 



In five dimensions, an orbifold was used to construct a supersymmetric model of (calculable!) 
electroweak symmetry breaking ^2j, predicting a light Higgs (related models were considered 



before in Refs. f2"3|, |25|). The realization of supersymmetry in this model is quite different from 
the usual minimal supersymmetric standard model — the theory below the compactification 
scale is not supersymmetric, while superpartners (the lightest being the stop) as well as mirror 
particles appear near that scale. In the six-dimensional case, a nonsupersymmetric T 2 /Z2 
orbifold construction was employed to build a higher- dimensional composite-Higgs model of 
electroweak breaking J23|. The issue of radius stabilization was not addressed in refs. [22], 
23] — the Higgs-dependent part of the Casimir energy played a crucial role in the analysis of 



electroweak symmetry breaking in [p^g] , but the radius was assumed to be fixed. 

Originally, orbifolds were invented in string theory [^, |27| as a tool for constructing new 
consistent string backgrounds. Modular invariance of the worldsheet orbifold CFT and/or 

, |30[ place severe constraints on their consistency as fundamental 



tadpole cancellation 



theories. The field theory orbifolds that we discuss in this paper cannot, at least at present, 
be derived from known string constructions. Nevertheless, their success with electroweak sym- 
metry breaking is appealing enough to warrant further study. We will only demand that the 



^We note that all field theory orbifolds discussed in this paper have some periodic fermions, hence the 
nonperturbative instability of the Kaluza-Klein vacuum towards decay into "nothing," pointed out in [^jj and 
recently discussed in p5[ is absent. 
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orbifold theories be consistent as low-energy effective theories (hence the requirement that the 
theory of the zero modes be anomaly free; for a recent discussion, see |[31|), valid up to some 
energy cutoff scale. Consistency of the field theory description demands then that the cutoff 
be (at least) an order of magnitude or so larger than the inverse compactification length scale. 

When performing loop calculations in the orbifold field theory, the orbifold boundary con- 
ditions lead to extra divergences at the orbifold fixed points, see, e.g., To cancel these, 
new terms localized at the fixed points have to be introduced in the lagrangian. These localized 
terms can be kinetic and mass terms, as well as interaction terms [|33], |33|]. Their coefficients 
are additional parameters of the orbifold theory. It is known ||35|| , that these "brane- localized" 
terms can affect the spectrum of the theory and consequently the Casimir energy. It is then 
natural to expect that they can play a role in the mechanism of radius stabilization. 

We will first work in an approximation neglecting the contributions of brane-localized kinetic 
and mass terms to the Casimir energy. This is a valid approximation if the coefficients of brane- 
localized kinetic and mass terms are of the order of the loop-generated values (consistency and 
naturalness demand that the tree-level coefficients be at least as large as the loop-induced values 
[|3G|1). We will then generalize the calculations to the case of larger brane-localized kinetic and 
mass terms. 

In Section O , we describe the five dimensional S l /Z 2 setup. We then present, in Section 
|2.2j , a general discussion of the divergences in the one-loop Casimir energy, the counterterms 
required for their cancellation, and the contributions to the radius potential from the Casimir 
energy and counterterms. The discussion of Section ^72] applies equally well to the 6-d case. 

Motivated by the 5-d and 6-d models, mentioned above, we calculate the Casimir energy 
in S 1 /Z 2 , in Section and T 2 /Z k compactifications, in Section f|. We compute the gravita- 
tional contribution, as well as those of even or odd massless matter fields obeying periodic or 



antiperiodic boundary conditions. In Section |27|, we also calculate the Casimir energy of fields 
with bulk mass terms in the 5-dimensional case and discuss a mechanism for stabilizing the 



radius with massive fields. The computations of Sections |2.3| , [2.4| , and f| are sufficiently general 
to be useful when considering radius stabilization in concrete five and six dimensional models. 

In Section |3], we include the contribution of the brane kinetic terms to the Casimir energy 
in the S 1 / Z 2 orbifold. We show that the brane kinetic terms, if they are sufficiently large, can 
lead to radius stabilization at a size bigger than the cutoff length scale. For example, a value 
consistent with "naive dimensional analysis" (NDA) can yield a radius several times larger 
than the cutoff length. NDA arguments |36j constrain the coefficients of brane-localized kinetic 
terms (with dimension of length) to be several times the cutoff length scale. 

A much larger value of the brane-localized kinetic terms can yield a radius much bigger than 
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the cutoff length scale. To investigate the viability of such a scenario as an effective low-energy 



theory, in Section |3.1| we ask whether the theory has a consistent perturbative expansion if the 
brane kinetic terms are significantly larger than the cutoff length scale. We investigate in detail 
the properties of the Green's function with the brane kinetic terms included and find that the 
brane-localized divergences are weaker, while bulk divergences are as in the theory without the 
brane kinetic terms; we note that this is consistent with the results of refs. |35], [H], [JSJ. The 



analysis of this Section strengthens the case for having values for the brane kinetic terms larger 
than those implied by NDA.Q 

Finally, in Section [5], we apply our results to the 5-d supersymmetric models of electroweak 
breaking |22j (see also f39j). We show that the Casimir energy contribution of the massless fields 



to the radius potential is repulsive (here we treat brane kinetic terms as small). We show that 
by fine-tuning the 5-d cosmological constant and "brane tension" counterterms (as discussed 



in Section |2.2|) , an acceptable minimum for the radius with vanishing cosmological constant 
can be achieved. The sign of the bulk cosmological term required to have such a minimum 
is positive, while the brane tensions need to be negative. It thus appears that the required 
bulk counterterm is not supersymmetric, at least in the simplest 5-d supergravities. This might 
present a problem, because a non-supersymmetric counterterm would introduce an additional 
source of supersymmetry breaking and potentially affect the predictions of the model. Thus, 
deciding whether this, or any other, stabilization mechanism is viable requires fully embedding 
the model in 5-d supergravity; this is a problem that we leave for future work. 

We present details of the calculation of one of the Casimir sums in an Appendix (other sums 
in the paper, where indicated, are computed similarly). We end with concluding remarks and 
a discussion of some open issues in Section 

2 Casimir energy on orbifolds. 

In this Section, we discuss the general issue of the Casimir energy on orbifolds, in the framework 
of the five dimensional S 1 jZi example. We begin by introducing the setup and notation in Sec- 



tion |2.1| . In Section |2.2| , we discuss the limits of applicability of the Casimir energy calculation, 
the divergences, and the counterterms required for their cancellation. The coefficients of these 
counterterms can not be computed from the low-energy effective theory alone — they are to be 
treated as parameters of the theory and are fixed by imposing normalization conditions on the 



# 2 This does not imply that the NDA arguments are incorrect. The main assumption of NDA — that there is a 
single fundamental length scale — does not hold when the brane kinetic terms take values larger than the cutoff 
length; we are merely saying that such a scenario can yield a consistent low-energy theory. 
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potential for the radius. Sections |2.3| , |2.4| contain the results of the calculation of the Casimir 
energy contribution of the gravitational field as well as of various massless and massive fields. 
In Section |2.4| we also discuss a possible stabilization mechanism using the Casimir energy of 
massive fields. 

2.1 The five dimensional S 1 fZi example. 

The general parameterization of the interval of an S 1 compactification is as follows, see, e.g., 

ds 2 = (g^ + A^A u (j)) dx^dx v + 20§ A^dx^dy + (f)^dy 2 . (1) 

Here the four- dimensional indices are denoted by /x and the five-dimensional coordinate y (taken 
to have dimension of length) is assumed to change between — L and L (until the physical radius 
is fixed, the scale L is a completely arbitrary length scale). For a general fluctuating background 
the fields g^, A^, and can depend on x v and y. The Z 2 orbifold is obtained after identifying 
points on the circle related by a reflection in the fifth coordinate, y ~ —y. The invariance of 
the interval (jl]) under the Z 2 symmetry determines the transformation properties of the fields: 
9nu(y) = 9v»{-y), A^y) = -A^-y), and <p(y) = <f)(-y). Since the field A^ is odd under the 
Z 2 symmetry, it can not have a zero mode. 

The parameterization of the metric in (||) is convenient, because in the four-dimensional 
effective field theory of the zero modes it gives rise to four-dimensional Einstein gravity with 
metric tensor g^ u (x), coupled to a dilatonp 3 ! 0(x) and, before the orbifold projection, an abelian 



gauge field A^(x). Since the Z 2 orbifold forbids the appearance of a zero mode of A^, we omit 
this field in what follows. The zero-mode effective theory is valid below the scale of the mass of 
the lowest Kaluza-Klein excitation, i.e. at energies below (0 1 / 3 L) _1 ; note that the physical size 
of the extra dimension is 1//3 L. More precisely, the five dimensional Einstein action, evaluated 
on a background ([!]) with the fields g^x), <f)(x), dependent only on x M is given by: 



Ml J d 5 xVGR{G) = Ml J d^xjg- 



R{g) + ^" 



6' 



(2) 



In the above formula G denotes the five dimensional metric tensor that can be read off eqn. ([!]), 
with A^ = and the rest of the fields only dependent on x M , while the four and five dimensional 
Planck scales are related by Mj = LM|; in the orbifold theory we only integrate over the 
fundamental region < y < 
# 3 We will occasionally call this field a 'radion.' 

^ 4 To avoid confusion, we note that the relation between four dimensional and five dimensional Planck scales 
involves the physical size of the orbifold, ^^L, instead of the arbitrary scale L; however, the form is more 
convenient before fixing (cf>). 
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2.2 Casimir energy, divergences, and count erterms. 

Our goal is to study the generation of a potential V(<f>) — and the possible existence of a 
minimum — for the dilaton field (j) due to quantum effects of the nonzero modes of the fields 
g,A,4>, as well as to quantum contributions to the Casimir energy of matter fields (massless or 
not) that might be present. 

We will perform the calculation of the Casimir energy around a constant background = 
Vnvi = 0, <f) = const., with 7/^-the Minkowski metric. Dynamical issues, such as the 
time evolution of the background as a backreaction to the Casimir energy, can then only be 
studied for time intervals such that the deviation of the metric from the assumed constant 
background is small. In this paper, we are interested in the existence of static stable minima 
with vanishing four-dimensional cosmological constant. Other interesting issues, such as the 
cosmological evolution of the background, are left for future work. 

We do not attempt to say anything about the cosmological constant problem here. To 
achieve vanishing cosmological constant at the minimum, we will resort to the fact that the 
computation of the Casimir energy is plagued with divergences, whose cancellation requires 
adding counterterms to the action (0). Divergences are short-distance phenomena and the 
counterterms needed to cancel them are local terms, respecting the short-distance symmetries 
of the theory. The Casimir energy, on the other hand, is a global effect, depending on the 
topology and boundary conditions of the compactification, and can not be described by a local 
term preserving the short-distance symmetries. For example, in the S 1 compactification, (bulk) 
counterterms (using a generally covariant regulator) should respect 5-d general covariance and 
should not depend on the fact that at large distances it is broken by the compactification. For 
the constant flat metric background of interest to us, there is only one divergent counterterm — 
the 5-d cosmological constant term (this is, strictly speaking, true only in the unorbifolded 
case, see discussion below). This term is: 



treat its coefficient a (of dimension (mass) 5 ) as a parameter to be fixed by the normalization 
conditions of the potential: 




(3) 



and thus contributes a potential ~ <f>' 



1 / 3 to the four- dimensional effective action (§). We will 



v{4>) = v'{<t>) = o . 



(4) 



We note that for values of 0, such that eqns. (H) hold, the metric background is flat and thus 
the calculation of the Casimir energy — which assumed that — is self-consistent. 
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It is clear that if only massless fields are present and the only scale in the problem is M5, 
the potential for due to the massless fields' Casimir energy is monotonic in and minimized 
either at = or 00, depending on the matter content of the theory. We will see that the 
massless fields' Casimir contribution to the effective action (0) is proportional to: 

-L / Af 2 (5) 

If the potential (|5]) is minimized for — > 0, one has to resort to strong curvatures or other 
effects beyond the reach of the effective theory to stabilize the radius and the effective theory 
description of the low-energy physics is not valid. If, on the other hand, the minimum is at — > 
00 one has a typical runaway potential (exponential, if the kinetic terms are rendered canonical 
by a field redefinition). This potential may provide for an acceptable form of quintessence. Our 
calculations show that in the five- dimensional models with the supersymmetric standard model 
in the bulk (with the minimal matter content) it is the second possibility that is realized. 

Adding the cosmo logical constant counterterm ~ -1 ^ 3 , eqn. @, changes the shape of the 
potential (|^) and can generate a minimum. However, it is impossible to have a vanishing cos- 
mological constant at the minimum with only the two contributions (|j) and (H) to the potential. 
The nonzero cosmo logical constant at the minimum renders the calculation inconsistent. 

Clearly, the resolution is to introduce another scale in the problem. In what follows we 
consider the possible ways to do this. 

The first, "nonminimal," way is to introduce additional fields in the theory and is applicable 
equally well to the unorbifolded theory. One can, for example, include a field with mass smaller 
than the cutoff of the five dimensional theory. In Section [2.4|, we calculate the Casimir energy 
of massive fields and show that upon adding the three contributions — massless, massive, and 
cosmological constant counterterm — it is possible to achieve a stable minimum for with 
vanishing cosmological constant, thus rendering the calculation self-consistent. The physical 
size of the extra dimension then is typically of order the inverse mass of the field. However, we 
will see that this mechanism, though in principle viable, can not be implemented in realistic 
models with the (supersymmetric) standard model fields in the five-dimensional bulk. 

A second possibility — we call it "minimal," since it does not involve introducing new fields 
solely for the purpose of radius stabilization — is to explore the fact that we are compactifying 
not on a smooth circle, but on an orbifold. The Green functions G(x — x';y,y') of fields with 
orbifold boundary conditions have, in addition to the usual short- distance singularity when 
x ti x ij,/^ y _^ yi ^ a singularity localized at the fixed points of the orbifold y = y' = 0, L (see for 



example ref. [[!2|], where the Green function of a scalar field in four dimensions with Dirichlet 
boundary conditions at a fixed two-plane is discussed). These singularities lead to divergent 
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terms localized at the fixed points; for recent discussions see refs. |p3| , p4 . 

For example, since the Casimir energy is the expectation value of the energy momentum 
tensor, related to the Green function (e.g., of a scalar field) by: 

(T MN (X)) = (d M d' N + . . .)G(X,X')\ X ^ X , , (6) 

(here X denotes both x^ and y), the singularities of G(X,X') at the orbifold fixed points lead 
to divergent contributions to the Casimir energy localized at the fixed points, in addition to 
the 5-d bulk cosmological constant term mentioned above. Canceling these extra divergences 
requires adding counterterms localized at the fixed points (i.e. "brane tensions"). We will treat 
the coefficients of these terms as parameters, which can not be determined in the low-energy 
theory. (We have to assume that the "brane tensions" are small enough in order to ignore the 
warping they typically produce; we will comment on this in our discussion in Section |5|.) For 
the metric background ([[]), these localized terms on the S 1 jZi orbifold are of the form: 

P I d 5 x5{y) Jg~ = /3 fd 4 x </r 2/3 , (7) 



where g denotes the induced metric at y = 0, and a similar term at the y = L fixed point. 
Thus the localized terms scale differently from the terms due to the cosmological constant 
counterterm ([3]) and the massless fields' "bulk" Casimir energy (|5]). We will show that in some 
cases it is possible to achieve a (local) minimum of V(0) with vanishing cosmological constant. 
The physical size of the extra dimension is then set by the coefficients of the tension terms. 

The singularity of the Green function leads also to other divergences, for example to kinetic 
terms localized at the orbifold fixed points [Q; these have been considered recently from a 



somewhat different perspective in refs. |35|, |37|, |38[ . The localized kinetic terms have the form: 



<fx c 5{y) dfidPQ (8) 

and introduce a new length scale, c, in the problem. The localized kinetic terms can significantly 
change the spectrum of the nonzero modes |35| , and hence the 0-dependence of the Casimir 
energy. We calculate, in Section |3], the Casimir energy in the theory with the inclusion of 
such terms and find that it is possible to achieve a minimum for with vanishing cosmological 
constant; we also discuss localized mass terms. 

Finally, even though we will not pursue it here, we should mention the possibility to combine 
the "minimal" and "nonminimal" approaches, exploring the fact that the compact space is an 
orbifold and at the same time introducing new fields for the purpose of size stabilization. This is 
much in the spirit of mechanisms considered before in codimension one and two in refs. pi f|T . 
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These are classical mechanisms: one postulates that the fixed points are sources for some bulk 



field, e.g., a scalar field. In codimension two, see f41| , one introduces two kinds of massless 
fields, "dual" to each other, and finds that the energy of one type grows with size while that 
of the other decreases (to achieve that, one has to impose certain boundary conditions at the 
fixed points), leading to a minimum for the size. Other "nonminimal" mechanisms have been 
discussed in pE2| , pE5[] . We note that if one does not want to achieve a hierarchically large radius, 



it is possible to use a classical mechanism of the type discussed above to balance a repulsive 
Casimir force. 

2.3 Gravity and massless matter fields. 

In this Section, we calculate the Casimir energy due to the gravitational field and to other 
massless fields, with both periodic and antiperiodic boundary conditions on the S 1 . Here we 
assume that the contributions of the various terms at the fixed points can be neglected, while 
Section ||| contains a discussion of the more general case. 

Most of the discussion in this Section is not new. For completeness, we briefly review the 



calculation of ref. |L3[ of the contribution to the Casimir energy of the fields g^ u , A^, 0. We then 
generalize the calculation to include matter fields with both periodic and antiperiodic boundary 
conditions on the S . Finally, we give the generalization to the orbifold case. We will see that, 
if the contributions from fixed-point-localized terms are negligible, the Casimir energy of the 
S l /Z 2 orbifold is one half that of the S l compactification. 

Later on, in Section [|, we apply the formulae obtained for various fields to calculate the 



potential for generated in the models of ref. p2|, where the supersymmetric standard model 



lives in the five dimensional bulk. We show that in these models the Casimir potential for 
is repulsive, hence (for vanishing coefficients of the counterterms) the compact space tends to 
expand to infinite size. 

We begin with the Casimir energy of the gravity sector. As in ref. [|T^], only the nonzero 
modes contribute to the potential for 0, and the contribution of the gravitational multiplet 
equals that of five massless real scalar fields. This can be seen after appropriate gauge fixing 
and is essentially due to the fact that the five- dimensional graviton has five polarization states 
PI - Thus, it is enough to compute the vacuum energy of a single massless scalar field, with 
periodic boundary conditions along the S 1 : 



2 n ±^ (2^) 4 "V <PL\ 
d 
ds 



= -4-C + ' SCalar ^)\s=o ■ (9) 



S 



The second line indicates that we use ^-function regularization to calculate the Casimir energy. 
Infinite contributions to the Casimir energy are thrown out by the regularization and countert- 
erms have to be added by hand, as discussed above. The periodic scalar (^-function is defined 
as: 

1 1 n 4 - 2s 

= 16vr 2 (2-s)(1-s) L 4 " 2 ^ 2 - C(2S ~ 4) ' (10) 

with — the Riemann zeta function. Plugging into eqn. (§) yields for the massless periodic 
scalar contribution to the potential in eqn. (Q): 

V+ ,s C alar (6) = ± = _ ^1 JL (U ) 

K9) ' 16 2 L 4 64vr 2 2 L 4 ' 1 ' 

which is the result of Jl3|| . We thus find that the contribution of the gravitational fluctuations 
to the Casimir energy, equal to 5V + ' scalar (<f)), is attractive (C(5) — 1.034), i.e. the circle tends 
to shrink to zero size. 

Having in our disposal the result for the periodic massless scalar, we can easily enumerate 
the Casimir contributions of all other periodic massless fields, e.g. by using knowledge about 
five-dimensional supersymmetry multiplets. Below, we summarize the results for all massless 
periodic fields of interest: 

y+,graviton^^ _ g y+, scalar fj\ 
y+Jermion^ _ _^ y+,scalar /^\ 

y+,vector^ = 3 y+,scalar^ ^ ^ 
y+,gravitino^±\ _ g y+,scalar r±\ 

It is also of interest to compute the contribution of matter fields (hypermultiplets) with 
antiperiodic boundary conditions on the S 1 . The scalar contribution V~' scalar ((f)) is given by 
the same expression as ((^) except for the sum being over half-integers. Thus, the formula for 
the antiperiodic zeta function ([H]) remains the same safe for the replacement ((2s — 4) — > 
((2s — 4, 1/2), where ((s, 1/2) = (2 s — l)C(s) is the generalized zeta function. Evaluating the 
derivative, we obtain a repulsive potential for <p from an antiperiodic real scalar field: 

y- scalar^ = y+, scalar ,^ (13) 

16 

while the contribution of an antiperiodic fermion is attractive: 

y-Jermion^ = \_ y +, scalar ^ _ ^ 
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Finally, we come to the S l /Z2 orbifold Casimir energy. Neglecting the effects of terms 
localized at the fixed points, the orbifold amounts to simply throwing out all modes even under 
y — > — y (for odd fields) or odd (for even fields). Since only the non-zero modes contribute to 
the Casimir energy, it is easy to see that in each case this amounts to throwing out half the 
modes in the sum over n in (|9]). (We note that an additional factor of 1/2 will occur because 
in the orbifold we only integrate from y = to y = L. This, however, will be irrelevant for us.) 

As noted in the Introduction, including the contribution to the Casimir energy of the mass- 
less fields alone (without including brane tensions) it is not possible to obtain a minimum for (f) 
with vanishing cosmo logical constant. We will, in what follows, use the results from this Section 
to discuss ways of obtaining a minimum when other possible contributions are also included. 

2.4 Casimir energy and radius stabilization with massive matter 
fields. 

In this Section, we calculate the contribution of both periodic and antiperiodic massive fields 
to the Casimir energy. We begin with the periodic case, when the vacuum energy for a real 
scalar of mass \i is given by: 

<#,,) - \tj^ ( k2 + ^ + jm)= -s<r-WU • («) 

The massive periodic scalar ^-function is: 

scalar /N 1 ^ f ^ f , 2 fV V? V 



2 n ^J (2tt) 4 V 0L 2 0V3 

1 



32vr 2 (2-s)(l-s) L 4 " 2s 2 - S 

where F(s; 0, c) is the series: 



F(s-2;O,/zL03/7r) , (16) 



F(s;0,c) = J2 (n 2 + c 2 )~ S , (17) 

n=— oo 

and can be evaluated as described in the Appendix. The contribution of an antiperiodic scalar 
field can also be evaluated — the only change is that the sum in flP7| ) is over half integers. 

The results for the Casimir contribution V +,scalar (<f), fi), V~ ,scalar ((f), n) of periodic or an- 
tiperiodic massive scalars respectively — omitting both finite and infinite contributions to the 
bulk cosmological constant — can be written in a common form: 

V ± ' scalar ((f),ii) = ° ' 



64tt 2 L 4 

/ i\ 

-. . / _•> / ,.,.,T 



Li 5 ( ±e~ 2L ^ ) + 2L^ Li 4 ( ± e - 2L ^ ) + ^L 2 /i 2 0§ Li 3 ( ±e~ 2L ^ 
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where Li n (x) are the polylogarithm functions. Since Li 5 (l) = £(5) and Li 5 (— 1) = — 15C(5)/16, 
the above expression reduces to the massless formulae for the periodic (pTTj) and antiperiodic ( |i~3| ) 
case in the small mass or small radius (<p — > 0) limit. As in the massless case, the contribution 
of a massive fermion on the S l is equal to minus 4 times the real scalar contribution; the Z 2 



orbifold contribution is one-half the unorbifolded contribution in each case. 

The bulk mass /x introduces a new scale in the problem, leading one to expect that stabilizing 
the radius using the Casimir energy of massive fields should be possible in certain cases. To see 
what the necessary conditions are, consider the following toy model with a single scale /i. Let 
there be a periodic field (scalar or fermion; or a number of fields) of mass /i and a massless sector 
(consisting of gravity and possibly other massless fields). Introducing the variable x = 2//L0 1 / 3 , 
we can write the total potential for x in the following form: 

V(x) = (ax* + (3 (hh{e~ x ) + x Li 4 (e^) + j Ii 8 (e-)) + 7) = x~ & f{x) , (19) 

where the first term is the cosmological constant counterterm, while the second and third terms 
are the contributions to the Casimir energy of the masssive and massless fields, respectively. 
That the potential V(x) can have a stable minimum is easy to understand qualitatively. A 
stable minimum can occur onlyQ if the massive contribution dominates near the origin (and is 
repulsive), while the massless term (with opposite sign) takes over at larger distances, where the 
massive contribution is exponentially suppressed. The normalization conditions (f|) become, 
for the potential V(x) (|T9|): 

/'(a;)=0, /(ar) = 0, (20) 

and can be shown to have solutions provided the coefficients 7 and f3 have opposite signs, while 
a and (3 have the same sign; in addition, the solution of (EO) is a minimum if 7 < 0. One also 
requires that | 7 | < Li 5 (l)|/3| so that the massive repulsive contribution dominates near x = 
ensuring the existence of a minimum. 

The conclusion one can draw from this simple one-scale model is that it is possible to have 
a stable minimum with zero cosmological constant, provided: i.) the massless contribution is 
attractive, ii.) the massive contribution is repulsive and dominant at small radii, and, in.) the 
bulk cosmological constant is fine-tuned (and positive, corresponding to a > 0); as usual, there 
is also a runaway minimum at x — > 00 -P 7 ! Barring unnaturally large coefficients, the radius is 
# 5 The potential eqn.(p"8|) was independently derived in Ref. p4| . 

# 6 It is important to note that the term multiplying (3 is a monotonically decreasing function of x equal to 
Lis{l) at x = 0. 

# 7 A similar conclusion can be drawn if one considers instead an antiperiodic massive field — once again one 
requires that its contribution to the Casimir energy be repulsive (hence the massive field is an antiperiodic 
scalar) and dominant at short distances in order to achieve a minimum (in this case a > as well). 
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stabilized at a size of order the inverse mass of the field, which is the only length scale in our 
toy model. 

Using the results of this Section, it is straightforward to include fields of different masses, 
as well as brane tension counterterms in the discussion of radion stabilization. 



3 Casimir energy contribution of brane- localized kinetic 
terms. 

In this Section, we calculate the contribution to the Casimir energy of a massless scalar field 
while taking into account the presence of kinetic terms localized at the fixed points y = 0, L. 
These terms can significantly affect the Kaluza-Klein spectrum of the field and change the 
Casimir energy. Readers interested only in the answer can skip to the final result, eqn. (|35|), 
and the following discussion of the radius dependence of the Casimir energy with boundary 
terms. 

In the metric background ([!]) with g^ v = T]^, = 0, = const., we take the quadratic 
action of a real periodic scalar field \l/ to bef^ 5 



l/d^x J dy d^d^ + ^(d y V) 2 + (2c <T 1/3 <%) + 2c L 0~ 1/3 5(L - y)) dfl&*V 



(21) 



where Co and cl are the two length scales introduced by the localized kinetic terms. Expanding 
the field ty(x,y) into Fourier components ^k(y) of four-dimensional mass y/Jfi, we find the 
equation of motion: 

d 2 + <\>k 2 + (2c o 0- 1/3 %) + 2c L <j)~ 1/3 6(L - y))# 2 ] * fc (y) = . (22) 

Now we introduce the variables k 2 = <pk 2 and Co,l = 4>~ l °o,Li and solve eqn. (B^) for even 
periodic fields *&(y + 2L) = ^(y), ^f(-y) = *&(y)- We find that the values of k for which the 
solution has the appropriate discontinuities to match the delta functions are determined by the 
solutions of: 

tan(jfeL) = - *fr + 5 *) . (23) 
1 - c c L k 2 

When co,l — ► 0, the equation simply gives the masses of the Kaluza-Klein modes \k n \ = \n\n/L. 

It is clear that the solutions of eqn. (p3|) are also labeled by an integer, and we shall denote them 

by k n . Note also that if c^l are positive, there are no negative k 2 solutions, i.e. no tachyons. 

# 8 Derivatives w.r.t. y will also generally appear. These require a more careful treatment in the thin wall 
limit we are studying here. For simplicity, we consider only the terms of eqn. (|2l|). These should be sufficient 
to illustrate the effects of localized terms on the Casimir energy. 
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Expanding the field ^(x, y) = { Pni.x)'^k n {y) in terms of the (even) modes ^f kn (y), we find 

n=0 

that the quadratic part of the four dimensional theory of the "Kaluza-Klein" modes (p n (x) is 
governed by the action: 

1 r 00 / k 2 \ 

- / d A X J2 d^nd^Vm + "f ¥>n¥>m Iran ■ (24) 

The normalization coefficients obey: 

L 

J dy^ kn (y)^k m (y) (1 + 2c %) + 2c L S(y - L)) = for m ± n . (25) 



L 

'Jmn 



-L 



The orthogonality of the Kaluza-Klein wavefunctions follows from the "Schrodinger" equation 



( p3| ) . The diagonal normalization coefficients 7 nn depend on n, however, in contrast to conven- 
tional Kaluza-Klein compactifications; this is important when considering interactions involving 
the Kaluza-Klein modes [[37], The normalization coefficients drop out of the Casimir energy, 
however. 

The Casimir energy of the real scalar field, even under the orbifold projection, can be 
then written as follows: 

1 ^_ r d 4 r> / _ P \ rl 

(26) 



^ C ^ c -)^E/0iog(^ + f) , - -c 

The ^-function is that of the operator in (|2"2"D and is given by: 



cq,l\ S J 



111 00 1 

( (s) = F(2s - 4) , where F(s) = Y — , (27) 

Sc °' lV ' 32vr 2 (2 - s)(l - s) L A ~ 2s (f) 2 - s y !l v ; ^x s n ' V 7 



where the sum is over the nonnegative roots of eqn. (p3|), which written in terms of x = kL, 
has the form taur = — xa/(l — bx 2 ). The coefficients a, b are expressed in terms of the ratios 
of the length scales co,l from ( |2"T| ) to the physical size of the orbifold l p h ys = Lcp 1 ^ 3 as follows: 

c + c L c + c L 
a = — - — = — , (28) 

^ ''phys 
C C L C C L 

b = ~1X = 1 2 ~- 

^ ''phys 

The sum F(s) of eqn. (|27|) can be taken by means of contour integrals as follows. Introduce 
the function: 

f(z) = az + (l- bz 2 )teaiz , (29) 
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and consider the contour integral: 

= I 'dz-J-^}- , C = ^ + C + + C e + C_ , (30) 
2m J z s j{zj 

where the contour is an infinite semicircle in the Kez > half-plane, C + runs along the 
imaginary axis from z = ioo to z = ie, C- — from z = —ie to —ioo, and C e is a small semicircle 
from z = ie to z = —ie in the Kez > half-plane. The sum F(s) in eqn. ([27]) can be written 
as: 

F (s) = Ei = !(*) + - s E 7 — = /(a) + ^— -CW , (31) 



n>0 X n Vr~ („.^i| 7T* 



where we used the fact that I(s) is determined by the poles of /'// inside the contour C, which 
occur at the zeros (z = x n ) and poles (z — (n+ 1/2)tt) of f(z), eqn. (|29|). 

We are interested in computing F'{— 4), which, as usual is defined by analytic continuation. 
The integral over vanishes for a sufficiently large positive Re(s). The rest of the integral 
can be explicitly evaluated and written as follows: 

s . (ks\°? _ s _ l ay + (1 + by 2 )tanhy 
I( s ) = - sin — / dy y log 



7r \ 2 / J '" 1 + ay + by 2 

e 

oo 

+ ^ sin ) J dy y^ 1 log(l + ay + by 2 ) (32) 

e 

jf ^e~ ise log/(ee ie ) - ^cos(y) . 



2tt 



-tt/2 



It is easy to see from (|32|) that the divergences as e — > cancel. This form of the integral is 
also appropriate for analytic continuation to negative s. The first integral is well-behaved for 
s < 0. The second integral can be explicitly evaluated in terms of hypergeometric functions; its 
analytic continuation for negative s vanishes as e — * 0. The analytic continuation of the last line 
for s < also vanishes in this limit. Thus, for negative s, we only have the first integral, where 
the limit e — > can be taken safely and we arrive at our formula for the analytic continuation 
of F(s) to negative s: 

t?( \ 28 ~ l t( \ , s ■ f ns \ 1a -s-u ay + (1 + 6y 2 )tanhy 

F \s) = — ({s) + - sin [—) dyy log — — — , 33 

tc s 7r\2/j 1 + ay + by z 

which gives rise to: 

15 7r 4 . r 3 ay + (1 + by 2 )t&iahy 



-£f(2.-4) 



s=0 
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In the limit a = b = 0, the expression (0), upon substitution in (p6|) , yields the Casimir energy 



lTD without the boundary terms. 

The final result for the Casimir energy of the even periodic scalar with boundary kinetic 
terms is, therefore: 



64vr 2 L 4 2 



L 0l/3 ' £202/3 j 



(35) 



oo 

with p(a,b) =4/ dyy 3 log [(ay + (1 + by 2 )tajahy) / (I + ay + by 2 )]. While the form of the <fi- 
o 

dependence is rather complicated and can be studied in detail only numerically, eqn. (KH) can 



be used to illustrate some qualitative features of the Casimir energy with boundary terms. 

Consider first the large-radius behavior — the limit — > oo. In this limit, the coefficients 
a, b vanish and eqn. (|35|) reduces to ( 11) — except for a factor of 1/2 due to the orbifold — giving 
thus an attractive potential at large 0: 

~ - + (*> 

Thus, as expected, the boundary terms do not affect the large-radius behavior of the Casimir 
energy. 

The limit of small radii <ft —>■ is more involved. To study it, note that the function p of 
can be written as: 



40 4 / 3 / dyyHog 



y + (1 + 6y 2 )tanh0 1 / 3 |/ 



1 + y + by 2 



(37) 



where the hat indicates that the <ft dependence has been scaled out of 6; in addition, to simplify 
our formulae, we have chosen the arbitrary length scale L such that a = 1. 

In the limit when the boundary term at one of the fixed points vanishes, i.e. b = 0, it is 
easy to see that p —>■ as <fi —>■ (the integral in (|37D diverges as <fi — * 0, but more slowly 
than the 4//3 prefactor, such that the whole expression (|37[) vanishes). The small-0 potential 
is therefore repulsive. Thus, in the small radius limit the Casimir energy for the periodic scalar 
with boundary term at only one of the fixed points is repulsive, while the large-radius behavior, 
eqn. (|36|), is attractive. We conclude that in the limiting case when the boundary kinetic terms 
at one of the fixed points, say cl = 0, vanishes, the Casimir energy exhibits a minimum for <fi. 
The minimum occurs for values of the physical radius of order the length scale cq set by the 
nonvanishing boundary kinetic term. 

The behavior of the Casimir energy with one nonzero boundary term can be inferred already 
from the eigenvalue equation (|23| ) in the limit when one of the boundary terms vanishes — then 
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the lowest eigenvalues are approximately the periodic Kaluza-Klein modes nir/L, while the large 
eigenvalues are the antiperiodic ones, i.e. (n + l/2)7r/L. One expects that the large-radius limit 
of the Casimir energy is dominated by the lowest eigenvalues — hence the attractive behavior 
characteristic of a periodic scalar — while the small-radius limit depends on the large eigenvalues, 
leading one to expect a repulsive behavior, as for an antiperiodic scalar. The behavior of the 
Casimir energy ( |2"6"D in the two limits <p — > oo and <fi — > confirms this expectation. 

In the case where both boundary terms are present, the analysis requires more care. Upon 
analyzing the integral (^) one finds that when a 7^ 0, b 7^ the Casimir energy is attractive 
at small as well (as in the case without boundary terms). Thus, the existence or not of a 
minimum depends on relative strength of a and b. By comparing to the a 7^ 0, b = case, one 
expects that for a sufficiently small ratio of b/a, a local minimum will still persist. A numerical 
analysis confirms this expectation — it is sufficient to have Cq/cl ~ .3 in order for a (local) 
minimum to exist. There is no minimum if Co = Cl, however. 

Thus, we conclude that the radius can be stabilized by the boundary kinetic terms, provided 
there is some asymmetry between the fixed points.^] As usual, the cosmological constant at 
the minimum can be adjusted to zero by tuning the coefficients of the counterterms. Finally, 
while for a periodic scalar with asymmetric boundary kinetic terms the minimum is local, we 
expect that for a periodic fermion the boundary kinetic terms will yield a global minimum. 

To conclude this Section, we consider the effect of brane kinetic terms on fields that obey 



antiperiodic boundary conditions. Solving the scalar equation ( f22|) but imposing antiperiodic 
boundary conditions: 

* fe (y + 2L) = -* fc (j/) , (38) 

we find that the modes can be divided in two types: i) modes that are even about y = L with 
a spectrum determined by: 

tan(jfeL) = , (39) 
c L k 

and ii) modes that are odd about y = L with a spectrum determined by: 

tan(A-L) = -L . (40) 
c k 

The Casimir energy of such modes is calculated in exactly the same way as for eqns. ( PU|) 

and (P7|). In the present case, the appropriate analytical continuation, analogous to eqn. 

# 9 We note that the classical mechanisms for radius stabilization of |^0|, |4lJ , also rely on an assymetry between 
the fixed points ("branes"). 
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yields: 



d 
ds 



F{2s - 4) 



15 7T 4 



C\-A)+A dy yHog^ 



7/tanhy + 1 



(41) 



s=0 
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■iL- l y + 1 



The zeros of the numerator inside the logarithm are those of eqns. (|39|) and (f40|), while the 
first term subtracts the contributions from the poles of the tangent. Individually, each of these 



antiperiodic field. In particular, for large <p the potential is repulsive. On the other end, for 
— > the spectrum is given by kL = nir, i.e. as for a field with periodic boundary conditions, 
which produces an attractive potential. Therefore, the potential must have a maximum at some 
finite 

3.1 Exact tree-level propagator including brane-localized terms 

In the previous Section, we showed that brane-localized terms can have an important effect 
on the Casimir energy and actually produce a stable minimum for the radion potential. This 
is not completely surprising since the coefficients of these terms introduce in general a new 
length scale in the problem. We showed that the inclusion of brane kinetic terms can stabilize 
the radion at a scale of the order of c^l of eqn. (|2l|) . Since these coefficients have dimension 
of length, and since we would like to stabilize the radius at a value somewhat larger than 
the inverse cutoff of the theory, it follows that the coefficients of the brane-localized kinetic 
terms need to be larger than the fundamental length scale if they are to be relevant for radion 
stabilization. In particular, their effects should be treated exactly, as we did in the calculation 
of the Casimir energy.^] 

Then the question arises as to whether it is consistent to allow these terms to be large, while 
treating other effects perturbatively. In this Section, we derive the exact tree-level propagator 
including brane kinetic terms and argue that perturbation theory does not break down, even if 
the quadratic brane operators have anomalously large coefficients. 

To this end, consider a real scalar field in 5-d flat spacetime with action: 



# 10 This behavior can be confirmed numerically from cqn. (f4l|); note that the integrand in (|4l|), after conve- 
niently choosing L — Cj, depends only on (p. 

# n One way localized kinetic terms could arise in, e.g., a string construction would be through the expectation 
values of fields, confined to propagate to the fixed points of the orbifolds ("twisted sector" fields); a similar 
mechanism has been exploited in a field theory context in |37[] . 
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In this Section, M = 0, 1,2,3, y, and /x = 0,1,2,3. We assume that the fifth dimension is 
compactified on S 1 and, to simplify notation, neglect the dependence on 0, which can trivially 
be included. The propagator obeys: 

(d 2 y + p 2 + 2c p 2 5(y) + 2c L p 2 6(y - L)) G(p; y, y') = 5(y - y') (43) 

where we Fourier transformed in the four noncompact coordinates. Since G(p; y, y') is periodic 
in y with period 2L, we can write: 

G( P ;y,y') = ^E^„.(m'), (44) 



s(y-y') = 



Replacing these expansions back in eqn. (f|^), dividing by p 2 — {^jf} and summing J2n e l ™ y , 
we obtain: 

G(p; y, y') + 2c p 2 G(p; 0, y')B(p, y) + 2c lP 2 G(p; L, y')B(p, y - L) = B(p, y - y') (45) 
where we defined: 

1 e~ iS r y 

B(p,y) = vr i: (46) 

2L n p 2 _ ™ 

Evaluating eqn. ( 45|) at y = and y = L gives two equations, from which G(p;0,y') and 



G(p; L, y') can be found as a function of B(p, y). Eqn. (|5|) then gives G(p; y, y'). So, it remains 
only to evaluate the sum in (0). This can be done by methods similar to those used to evaluate 
the sums for the Casimir energy (see Appendix) and we shall not repeat them here. The result 
is: 

%L) , (47) 



B(p, y) = j- csc(pL) cos ^1 - j. + 2 



y 



2L 

where [x] denotes the integer closest to x, but smaller than x. 
The propagator can then be written as follows: 

G(p; y, y') = G + (p; y, y') + G_ (p; y, y') (48) 

where G + (p;y,y') is the propagator for the even modes of the field and G_(p;y,y') is the 
propagator for the odd modes, that is (denoting $±(p, y) = ($(p,y) ± —y))/2): 

G±(p;y, y') = ($±(p,y)®±(-P,y')) (49) 

= 7 (G(p; y, y') ± G(p; -y, y') ± G(p; y, -y') + G(p; y, ?/)) • 
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The final expression for the propagators of the even and odd modes is easily obtained from 
eqns. ([45]) and (f47|). In the region y,y' G [0, L], the propagator of the even modes of the field 
is: 

G {1 \p-yy') = - [ cos (^<) ~ c oP sin (py<)} [cos(p(L - ?/>)) - c L psin(p(L - y>))] 
+ ' ' 2p [(c c L p 2 - 1) sin(pL) - (c + c L )pcos(pL)] 

where y K (y>) denotes the smaller (larger) of y, y' and p = \fjP. We observe that G + (p; y,y') 
has poles whenever (cqClP 2 — l)sin(pL) — (co + cl)p cos(pL) = 0. This agrees with our result 
from Section [3|, see eqn. (p3|). In the same region y' G [0,L], the propagator for the odd 
modes has the much simpler form: 

G { l\p;y,y') = csc(pL) sin(py<) sin(p(L -y>)) . (51) 

We note that G-(p;y,y') has poles at p = mr/L, which is consistent with the fact that the 
odd modes vanish at the location of the branes and therefore do not feel the localized kinetic 
terms. The orbifold projection amounts to keeping only G + (p; 0,0) for even parity fields and 
G-(p; 0, 0) for odd parity fields. 

The propagators in the other regions can be easily found from the reflection properties of 
$±(p, y) in y according to eqn. (|49|). Explicitly, when y, y' G [— L, 0], 

G ( i\p;y,y')=G { l\p;-y,-y'), (52) 

while for y G [0, L] and y' G [-L, 0] or v. v., 

Gf{p-y,y') = ±G { l\p-\y\Ay'\) ■ (53) 

For y, y' outside the [— L, L] interval, G±(p; y, y') are extended periodically in both y and y' . 

We also note that G+\p; y, y') — > G ( l\p; y, y') as both c , cl — > oo. It is easy to check, that 
in this limit the full propagator G(p;y,y') is precisely the propagator with Dirichlet boundary 
conditions. This can also be understood by considering the eigenmodes of the scalar field 
equation and noting that both the even and odd modes vanish at the branes when cq, cl — > oo. 
Thus, the limit of very large kinetic terms reduces to a problem where the branes act as perfect 
mirrors. 

We are now ready to address the main goal of this Section — analyzing the divergences 
occurring in the perturbative expansion when the quadratic brane-localized terms are large, 
while interaction terms are considered as perturbations. The expressions (BO) and fl5"T|) allow us 
to study the behavior of the propagator as a function of the 4-d momentum p and thus analyze 
the divergences that will appear in loops with this propagator. 
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We find, that if y ^ y' the propagator vanishes exponentially for large Euclidean momenta, 
p ^> L^ 1 . This shows that the possible divergences in the theory are local in y, as expected. 
When y — y' there is a difference whether we are sitting inside the bulk or at one of the branes. 
In the former case, G(p;y,y) ~ ^- which corresponds to a 5-d behavior. In the latter case, 
when y = y' at one of the branes, say y = 0, G(p; 0, 0) = G + (p; 0, 0) ~ 2 , + l Cop 2 with cl in place 
of Co when y — y' — L. This is precisely the behavior found in ref. |35| in the case of gravity.* 12 
On the other end, for p <C L" 1 , the propagator behaves like G(p; y, y') ~ ((2c + 2cl + L)p 2 ) -1 , 
which shows the 4-d behavior expected from the finite size of the extra dimension. 

If we add interactions to the theory, we can proceed to calculate Feynman diagrams with 
our exact tree- level propagator, eqn. (|48l), or the appropriate one in an orbifolded theory. We 
imagine that the theory is valid below some scale A, and that the 4-d momentum integrals are 
cutoff at A. We also assume that the radius has been stabilized at L ^> A -1 , so that there 
is a region of energies where the 5-d effective theory is valid. The asymptotic behavior of the 
propagator for large Euclidean momenta, discussed in the previous paragraph, then shows that, 
in the bulk, the possible divergences of Feynman diagrams are the same as in an uncompactified 
5-d theory, independent of the size of the brane kinetic terms. The convergence properties of 
Feynman diagrams at the branes, however, depend on the size of Co, Ct,. For definiteness, let us 
concentrate on the y = brane and write Cq = cA -1 . If the brane couplings are perturbative, 
i.e. c % 1, the high energy behavior of the propagator (~ i) is as expected in a 5-d theory. 
If, on the other hand c ^> 1, the convergence properties are improved — they are as in a 4-d 
theory — since for p ^> Cq 1 = A/c, the propagator behaves like ~ jp. Thus, it is consistent to 
assume that q ^> A -1 , in which case the brane kinetic terms can play a role in the stabilization 
mechanism for the radion. If one is not interested in obtaining a large hierarchy, moderately 
large q might be enough. We remark that, as shown in ref. the presence of brane kinetic 
terms seems to be quite generic in theories with orbifold fixed points. 

To conclude this Section, we note that it is straightforward to include bulk and brane mass 
terms as well. Consider adding to the action of eqn. fl42|) the terms: 

AS = -^Jd A x J L dy [M 2 + 2m %) + 2m 1 5{y - L)j $ 2 . (54) 

The analogs of eqns. ( p0| ) and ( |5"ID for the propagator derived from S + AS are: 

(1) , _ [cos(py<) - bop' 1 sin(py < )] [cos(p(L - ?/>)) - b L p~ x sin(p(L - y>))] 
+ {P,V,y> 2p[(6 6 L p- 2 - 1) sin(pL) - (b + b L )p~^ cos(pL)] 



# 12 The models of pj| do not assume that the extra dimensions are compactified. Their point is that the above 
4-d behavior due to the brane kinetic terms can make the compactification unnecessary. The fact that our 
result reduces precisely to theirs in the large momentum region is clear, since we do not expect the large energy 
behavior to be sensitive to what may happen far away from the brane in question. 
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G W (p;y,y') = - — csc(pL) sin(py<) sin(p(L - y>)) (55) 
ip 

where now hi = Cip 2 — m, and p = \/p 2 — M 2 . These terms also induce a nontrivial radion 
potential. 

4 Casimir energies in six dimensions compactified on a 
torus 

Now we consider the issue of the Casimir energy and radius stabilization by quantum effects in 
the richer case of a 6-d spacetime. We consider compactifying two of the dimensions on a torus, 
where the calculations run in parallel with our 5-d analysis. The toroidal compactification is 
the one employed in the composite-Higgs models of ref . . The detailed spectrum and current 
experimental constraints on the size of the T 2 /Z 2 orbifold from precision electroweak observables 



are discussed in ref. Additional motivation for considering orbifold compactifications of 

this type is the observation [^] that in 6-d standard-model global anomaly cancellation places 
severe constraints on the number of identical generations, requiring n g = mod 3. 

The main difference from the 5-d Casimir energy is that there are now three moduli, which 
would have to be stabilized. As we will see, the dependence of the Casimir energies on these 
moduli is nontrivial. To define a torus, one specifies a lattice in the plane and identifies points 
that differ by a lattice vector. The three moduli can be identified with the lengths of the two 
lattice vectors and the angle between them. Alternatively, we can use the area of the torus and 
the modular parameter r = n + IT2. We consider the background interval: 

ds 2 = A~ x ri^dx^dx" + A i ij dy i dy j (56) 

where y % G [0, L] and the metric 7^ on the torus is: 

•* = h( \ M> ) ' (57 > 

In this parameterization, det(j) = 1, so A is indeed the (dimensionless) area of the torus. 
Further, the parameterization is such that for fluctuations of the metric with 77^ — > g^^x^), 
A(x^), r(x /i ) depending only on the four noncompact coordinates, we have: 



Ml J d 6 xVGR{G) =M 2 J d A x^g~ 



R{9) + ~^r~ + a 2 



(5? 
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where M 2 = L 2 M,f p T3 | , R is the Ricci scalar evaluated with the metric indicated in parenthesis 
and f = T\ — 1T2. There are also two massless spin-one fields with their corresponding Kaluza- 
Klein towers. These modes also contribute to the total Casimir energy, but they do not play 
a role in the calculation of the Casimir energy itself (see below) so we have not written them 
explicitly. We will see below that their zero modes are projected out in the orbifold. 

As in the 5-d case, it is enough to calculate the Casimir energy of a real scalar field in 
the desired background (see discussion of gauge fixing below). Thus, we begin with the vac- 
uum energy of a single massless scalar field propagating in the background Q5BD, with periodic 
boundary conditions: 

= -^( + ' scalar (s)\ s=0 , (59) 
as 

where R = L/(2ir). Here we used again ^-function regularization to define this expression, 
where now: 

C+ ,scalar,, = 1 y / f , 22 2 , | 

{ J 2^ n J (2tt) 4 V 2 1 1 

7A(s-2)(s-l)^ n '\m-nr\^ ' (6 ° } 

and the prime indicates that the zero mode is to be excluded. It is clear from the second line 
that £+> scaiar ( s ) i s modular invariant, reflecting the fact that the Casimir energy is independent 
of the discrete choice of lattice vectors defining the torus (these choices are related by SL(2, Z) 
transformations of the modular parameter r, see, e.g., |^6|). We defer the evaluation of the 
sum in (|60|) to the Appendix. The result for the potential is: 

V + ' scalar = -^- 4 A- 4 f(r,f) , (61) 
with (taking T2 > without loss of generality): 



f(rf) - ^L r 3 + 3C(5)1 + 



4 00 / 31 31 \ 

-2 E K Li 3(<? m ) + tt- - m Li 4 (g m ) + — -j Li 5 (g m ) + h.c. , (62) 



# 13 M is not yet the physical 4-d Planck mass that characterizes the gravitational interactions. Once A gets 
a vacuum expectation value, it will be convenient to do a constant Weyl rescaling (j^u = A~ l (jfiv It is g that 
couples to all fields. As a result, it is the coordinate independent combination AL 2 that appears everywhere 
and the physical 4-d Planck mass is Mj = (AL 2 )Mq. 
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where Li n (x) are the polylogarithm functions and q = e 27r4 ( r i+* r 2). 

A few comments are due on the properties of the function /(r, r) in (|61|). As mentioned 
above, /(r, f) is invariant under SL(2, Z) modular transformations (r — > (ar+6) / (cr+d), where 
a, 6, c, d are integers, obeying ad — be = 1) due to the fact that tori with modular parameters 
related by an SL(2, Z) transformation are diffeomorphic. Distinct tori have modular parameters 
r taking values in the fundamental region |t| > 1, 1/2 > T\ > —1/2, r 2 > 0. It is therefore 
enough to study the behavior of the Casimir energy floT| ) in the fundamental region. The 
points r = i and r = e 2m ^ 3 of the fundamental region are fixed points of the transformations 
t — > — 1/t and r — > — 1/(1 + r), respectively, and are thus extrema of /(r, r). Numerically 
we find that the point t = i, corresponding to a rectangular torus, is a saddle point, the point 
r = e 2?ri / 3 , corresponding to a torus of angle 120 degrees, is a minimum, and that there are no 
other extrema of /(r, r). The function /(r, r) is positive definite and unbounded as r 2 — > oo 
in the fundamental region. Therefore, regarding the dependence of the Casimir energy on the 
modular parameter r, for fixed area A, we can conclude the following: i.) if the Casimir energy 
is attractive, as for a periodic scalar, eqn. (|6lD , the energy is unbounded below as r 2 — ► oo. 
ii.) if the Casimir energy is repulsive, the torus with r = e 2m ^ 3 minimizes the potential for a 
fixed area. For repulsive total Casimir energy, therefore, stabilizing the area modulus leads to 
stabilization of the modular parameter. 

To calculate now the total contribution to the Casimir energy due to the gravitational 
fluctuations, it is necessary to count the number of physical degrees of freedom. We write the 
6-d metric as Gmn = Gmn^^mn, where G^) N is the background defined in eqn. ( [56] ) and Hmn 
are the quantum fluctuations. To begin, note that in a flat background it is always possible to 
fix a "transverse-traceless" gauge: 

G^ MK d K h MN = 
h = 

where h MN = Hmn — \Gmn^ 1 an d h — G^ mn Jimn- This gauge fixing conditions elim- 
inate 7 of the original 21 polarizations. The fluctuations then satisfy the wave equation 
G^ KL dKdLh,MN = 0. However, there is still a residual gauge freedom parameterized by 
h' MN = h MN + d M £,N + d N £ M , where G ( -°^ KL d K d L ^ M = 0. Working in Fourier space (dis- 
crete along the two compactified dimensions), it is easy to see that, for modes with a nonzero 
momentum in the compact dimensions, there is enough freedom to set: 

h^ix^y 1 ,y 2 ) = -h^(x } y l ,y 2 ) 
h5(x,y l ,y 2 ) = h^(x } y l ,y 2 ) . 
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So the massive modes (from a 4-d perspective) include a tower of spin-2 fields (with 5 polar- 
izations each), a tower of spin-1 fields (with 3 polarizations each) and a tower of real scalars. 
All of these give a contribution to the Casimir energy which is 9 times that of eqn. (|61|) . 

The zero modes include (before orbifolding) the 4-d graviton, two gauge fields, and the three 
moduli that define the size and shape of the torus; as usual, the zero modes do not contribute 
to the Casimir energy. We also quote the contributions of other relevant bulk fields: 



y+j graviton _ g y+i scalar 
+ ,W eyl j 'ermion ^ y 



) 

scalar 



y+,vector ^ y+i scalar (63) 

y+,2form g y+, scalar 

y+,gravitino y+, scalar 

The simplest way to arrive at these results is by knowing the 6-d (1, 0) supersymmetry multi- 
plets. The toroidal compactification preserves supersymmetry, hence the Casimir energy should 
vanish within supermultiplets. The second equation follows from the fact that the (1,0) hy- 
permultiplet contains a complex Weyl fermion and two complex scalars, while the third follows 
from the content of a (1, 0) vector multiplet (a vector field and a Weyl fermion). The self-dual 
two form field contribution is denoted by y+> 2 f° rm (recall that a (1, 0) tensor multiplet consists 
of a Weyl fermion, an anti-self-dual two form antisymmetric tensor field, and a real scalar). Fi- 
nally, the last line follows from the structure of the (1,0) gravity multiplet: graviton, gravitino, 
and a self-dual two form. 

Note that if we orbifold the torus by a discrete subgroup of the 2-d rotation group, it 
is possible to "freeze" the values of T\ and Ti and concentrate on the A dependence of 
alone. ^ For example, if we write the coordinates in the extra dimensions in complex notation 
as z = y 1 + iy 2 and identify points z ~ e 2nt ^ n z, it is easy to see that whenever n > 2, the 
requirement that the lattice that defines the torus have the appropriate symmetry implies that 
T\ and T2 must take on particular values.^ 

The orbifolding will also project out some of the modes of the gravity multiplet as can 
be seen by the requirement that the line element be invariant. For example, the case n = 4 
corresponds to the identifications (y 1 -,y z ) ~ (— y 2 ,y x \ Writing the line element as: 

ds 2 = G^ u dx^dx u + A fli dx >l dy i + 'j ij dy i dy j , (64) 



^ 14 This would probably be necessary in any case to obtain chiral fermions in the low energy theory. 

# 15 In particular, the case n — 4 corresponds to a square torus with t\ = 0, T2 = 1. In this case, the sum 
^2'(n 2 +m 2 )~ s in eqn. (|60|) can be evaluated more easily with the help of the Jacobi identities. We have checked 
that eqn. (|6^) reproduces this result. 



24 



determines the transformation properties of the various metric components: 



G,„(x»,-y 2 ,y l ) = G, u {x^y\y 2 ) 

A^x^-y^y 1 ) = -A, 2 (x^y\y 2 ) 

A^{x»,-y 2 ,y l ) = A^{x^y\y 2 ) (65) 

7n(a; M ,-?/ 2 ,y 1 ) = 722(3:", y\y 2 ) 

712(2^, -y 2 , j/ 1 ) = —721(2^, y 1 , y 2 ) 

These parity assignments imply A^x**, — y 1 , — y 2 ) = —A^ix^^y 1 ^ 2 ) and 7i2(a: / \ — y 1 , — y 2 ) = 
— 712(2^, y 1 , y 2 )- So the zero modes of these components are projected out. Also all their 
Kaluza-Klein components are forced to vanish at the positions of the orbifold fixed points 
(y 1 ,?/ 2 ) = (0, 0), (0, L), (L, 0) and (L,L). In addition, the other (y-independent) zero modes 
(apart from the 4-d graviton) satisfy 7n(x /1 ) = 722(2^), so there is only one modulus left, which 
can be identified with the area of the torus. 

Quite generally, there will be a discrete set of points that are left invariant by the orbifold 
identification. As discussed in the previous sections, we expect the presence of terms in the 
action localized on these fixed points. Here we consider only the effect of the brane tensions, 
which gives a contribution to the potential for the area modulus A: 



T J d G x 5{y) J Z g = T j d A x A' 2 (66) 

where g is the induced metric at the relevant fixed point and T — the corresponding brane 
tension. Adding these terms, as well as a bulk cosmo logical constant: 

A / d 6 x VG = AL 2 f d 4 x A' 1 , (67) 



to the contributions of the massless fields (|BT ) it is possible to obtain a minimum for A that 



satisfies the normalization conditions (01), as in the 5-d case. 



If we sum the contributions of all fields to the Casimir energy in the model of ref. [23| 



we find that the potential for the area modulus is repulsive. We note that ref. |23| used a Z 2 



orbifold, in which case there are two noncompact moduli; it is, however, possible to generalize 
the construction to higher Z^s, affecting only the massive spectrum, and thus project all but 
the area modulus (alternatively, since the Casimir energy in the model with minimal matter 
content is repulsive, as remarked after eqn. (|61|), for fixed area the torus with r = e 2nl ^ 3 has 
minimal Casimir energy). With just the area modulus, a minimum with zero cosmological 
constant can be obtained after fine tuning the coefficients of the counterterms fl66|), flBTp. 



25 



5 Applications 



In this Section, we consider the question of radius stabilization in the supersymmetric 5-d model 
of [^2J (see also [p^fl ). This model is highly predictive, and the authors of ref. |22j were able to 
determine not only the Higgs and superpartner masses, but also the compactification radius of 
the fifth dimension, which, in fact, sets the scale for all these masses. However, they did not 
provide a stabilization mechanism for the radius. We are interested in determining whether the 
radius can be stabilized by quantum effects due to the original (minimal) field content. 

Let us briefly review the main elements of the model of ref. f22[. One starts from a five 
dimensional setup where the standard model superfields propagate in the bulk. The fifth 
dimension is compactified on S , which is then orbifolded by identifying points that are related 
by two Z 2 parities. The first one gives just the S 1 jZ 2 orbifold studied in Section || and breaks 
half the supersymmetry. The second identification breaks the rest of the supersymmetry giving 
precisely the standard model at low energies (below the compactification scale). 

The standard model gauge fields are described by 5-d vector supermultiplets (A M , A, A', a), 

which from a 4-d point of view decompose into one N = 1 vector supermultiplet V = (A 11 , A) 

and one chiral superfield in the adjoint representation E(0s,'0s), where $s = (cr + iA 5 )/\/2. 

The standard model matter and Higgs fields are described by hypermultiplets (\1/, $, $')x? where 

X = M (= Q,U, D, L, E) for matter fields and X = H for the single Higgs field of the model. 

Here \1/ is a Dirac fermion and $, are two complex scalars. In 4-d language, these correspond 

to two N = 1 chiral supermultiplets X($x,V ; x) an d X c ($ x ,ip x ), where \l/ = (ipx^x) an d 
$' = $ c t. 

In this model, matter and Higgs superfields are distinguished by different parity assignments 
under the Z 2 x Z ; 2 projection. An equivalent way to describe this projection, which is more 
suited to the notation employed in Section |2.3| is the following: we start from the S 1 /Z 2 orbifold 
discussed in Section |2| and break the remaining 4-d N = 1 supersymmetry by the Scherk- 
Schwarz mechanism, imposing 

$(y + L) = R P $(y), (68) 

where Rp is the R-parity of the unbroken supersymmetry of S 1 /Z 2 and $ represents any 
component field. Eqn. (|68| ) is the same as imposing periodic (antiperiodic) boundary conditions 
on the fundamental interval [0, L] on the fields with Rp even (odd). If we denote the possible 
parity assignments by (+, +), (+, — ), (— , +), (— , — ) where the first entry refers to the Z 2 and 
the second to the Rp parity, the component field parity assignments are as given in Table 1. 
In addition, one can write supersymmetric interactions at the orbifold fixed points (using 
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(+,+) 


( + ,") 


(-,+) 


(- -) 


M 










H 


<$>H 




<T) c t 




V 




X 







Table 1: Parity assignments for matter, Higgs and vector fields in the model of ref. [fffi| . The 
first assignment refers to the Z 2 parity and the second to the Rp parity of each component 
field. 



the notation of ref. f22| ) : 



AL 



~(8(y) + 8(y-2L)) 



d 2 6{XuQUH) 



L) + 5{y + L)) J d 2 6\X D Q'D'H' + X E L'E'H') + h.c. 



(69) 



which give rise to standard model quark masses after electroweak symmetry breaking. 

Ref. |22 calculated the one-loop effective potential for the standard model Higgs taking into 
account the contribution of the top Yukawa coupling (the other Yukawa couplings are much 
smaller and are neglected). We note that this is precisely the Casimir energy due to the Q 
and U matter superfields in the presence of brane masses (= Xu{H)). Adding the tree-level 
potential (from the D-terms) to this quantum contribution, the authors of ref. |22| found a 
Higgs potential of the form (see eqns. (15) and (31) of ref. [f22]]): 

1 



BHN \ 



H) 



l L 4 



V(<P 1/3 LH) 



(70) 



where is the radion field in the parameterization eqn. ([]]). The minimization of this potential 
with respect to H fixes the combination x = cj) l ^LH at a value xq that can be determined 
numerically. Knowledge of the Higgs vev then determines the size of the extra space 
However, eqn. ([Tf]) is only part of the contribution to the potential for 0. As discussed in Section 
0, all other fields, as well as possible cosmological and brane tensions contribute to V(<j)). 

To analyze this potential, we first neglect the possible contributions from brane kinetic 
terms. Then the spectrum for each of the parity assignments is: 



ni 



rn 



m 



(+,+) 

n 

(+,-) 
n 

(-,+) 



m 



-) 



2mr/L 
{2n + 1)tt/L 
(2n + 2)vr/L 
(2n + 1)tt/L 



(71) 



with n — 0, 1, 2, We can obtain the Casimir energies in the present orbifold from eqn. ( |T0D 

by making the replacement L — > L/2 and including a factor of \ due to the first Z 2 orbifold, 
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thus giving a total factor of 8. Thus, the contribution from the KK tower of each real degree 
of freedom is: 



y(+,+),realscalar / 
y (+,—), realscalar / 

y(-,+),realscalar(±\ _ gy +, scalar (±\ 
T/(— >— ),realscalar / 



(0) = 


8^ H 




(0) = 


8V~ 


-, scalar 


(0) = 


8V A 


-, scalar 


(0) = 


8V 


-, scalar 



where ^/+> sca/ar (0) and \Z-> sca/ar (0) are given in eqns. (11) and fli"3|), respectively. Taking into 
account the number of degrees of freedom and parity assignments given in Table 1, we get for 
matter, Higgs and vector supermultiplets the following contributions to the scalar potential: 

V {M \(t)) = -62V + ' scalar ((p) 

V {H \(j>) = +62V + > scalar (0) (73) 
V (v \0) = +62V + ' scalar (0) 

If we consider now the field content of the present model, namely the 577(3) x SU(2)l x U(1)y 
vector multiplets, three generations of matter hypermultiplets (less the third generation Q and 
U superfields which were included in eqn.([7D|)), and one Higgs hypermultiplet, we get a total 
contribution: 

rmassless/i\ 1 oe / i /•+, scalar i i 

. , 1023C(5) 1 
16 

Thus the total radion potential is: 



ymassiess,^ = -13Q4V + ' scalar = + — ■ (74) 



LV(0) = a<P~ 1/3 + (3<P~ 2/3 + 1<P' 2 , (75) 

where 7 = V(x) + 102 ^^ 5 ) includes the contributions from all matter fields and we also included 
a bulk cosmological constant a, and possible brane tensions (3. As we said before, the min- 
imization with respect to H, V'(x) = 0, fixes x = <p l ^LH at some Xq and also the value of 
V = V(x ) in eqn. (|70|) (which turns out to be positive). The equation of motion for then 
requires: 

- r 1 (f r 1/3 + f 0- 2/3 + 2 7 0- 2 ) = . (76) 

In addition, the requirement that the cosmological constant vanishes implies: 



ac 



We note that solving eqns.([HJ) and ( |77|) imposes a fine-tuning condition among a, [3 and 7. 
This is of course the cosmological constant problem, which we are not trying to address here. 



28 



Note also that we found that 7 > 0, which is a consequence of the fact that the matter content- 
which produces a repulsive potential — dominates over the gauge and Higgs contributions.!^ 



We find then, that eqns. (|76|) and (|77[) have a solution — which is a minimum — provided a > 
and (3 < (and the fine-tuning condition is satisfied). Thus, this mechanism could only 
stabilize the radius if the bulk cosmological constant was positive (de Sitter) and there are 
negative tension branes. The first condition seems to be especially problematic, since a positive 
cosmological constant seems to be incompatible with supersymmetry. We also note that the 
minimum so obtained is not the only minimum. Since the potential vanishes for large radius, 
there is a second degenerate minimum at <fi = 00. 

It is possible — and we leave this for future p 17 ! study — that the contribution from brane 
kinetic terms that we discussed in Section |3| could change the potential enough to produce 
a minimum with zero 4-d cosmological constant even when the bulk cosmological constant is 
negative. However, we note that this would be at best a local minimum, and the true minimum 
would be anti-de Sitter. To see this, we note that the cosmological constant term always 
dominates at large separation of the branes (i.e., at large 0), so for negative bulk cosmological 
constant the potential will always tend to zero from below. If there is a minimum at a finite 
<p where the potential vanishes, there will always be a second minimum where the potential is 
negative. In this case one would have to worry about the tunneling probability to the true AdS 
vacuum (recall, however, the suppression of tunneling to AdS space f47|j ). 

To conclude this Section, an issue which deserves a comment is that of the scale of 5-d 
gravity (and of 6-d gravity from the previous Section) in this model and the consistency of our 
approximations, which ignored warping due to brane tension and bulk cosmological constant 
counterterms. This is clearly a concern, since the 5-d gauge theory is non-asymptotically 
free and the couplings blow up a decade or so above the compactification radius, demanding 
a transition to a more fundamental description of the theory somewhere in the (multi)-TeV 
region. If the fundamental scale of 5-d (or 6-d) gravity was in that range as well, we would 
need a mechanism explaining the smallness of the observed gravitational interaction in four 
dimensions. 

One possibility in the present context is to assume that the higher dimensional gravity scale 

is in the multi-TeV range while the weakness of 4-d gravity is due to the presence of additional 

# 16 And over the gravity multiplet — it is easy to check that its attractive contribution, neglected in (^), is 
much smaller than the matter fields' repulsion. 

* 17 It would be also interesting to see to what extent the predictions of the model depend on the values of the 
brane kinetic terms. In particular, bounds on the coefficients of the kinetic terms from single Kaluza-Klein- 
mode production should be revisited — as follows from eqns. ( p4| , |25| ), the interactions of these modes depend 
nontrivially on their excitation number. 
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large dimensions (e.g., of millimeter size) accessible only to gravity, in the spirit of 0. The 
radius stabilization mechanism is then more complicated due to the presence of these extra 
dimensions (also, some of our formulae for the Casimir energies would need to be modified). 
More of a concern in this regard is the neglect of the backreaction of brane tensions and 
cosmological constant counterterms — since their size is determined essentially by the inverse 
compactification scale (of order TeV), which is not much smaller than the scale of the higher 
dimensional gravity. 

Another possibility — which makes neglecting the backreaction of branes and bulk more 
palatable — is to assume that the 5-d (or 6-d) gravity scale M is close to the four dimensional 
observed Mpi anc k and the two are related in the usual way, e.g., Mp lanck ~ L p h ys M 3 for the 
5-d case. From a low-energy point of view this constitutes an unexplained fine-tuning, but 
may have its origin, as recently pointed out in ||48|| , in "little string theory" (the term first 
appeared in [jl9 |; for a review, see J5(|). A major assumption required in order to embed the 



models discussed in this paper in this framework is that there is a window of energies above 
the compactification scale where the field theory description is still valid. 

6 Concluding remarks. 

This paper was devoted to the Casimir energy in five and six dimensional field theory orbifolds 
and its effect on radion stabilization. We were motivated by recent models of electroweak 
symmetry breaking, which used 5-d (supersymmetric) E2] and 6-d EM orbifolds. 



We gave a general discussion of the divergences of the Casimir energy and the countert- 
erms required for their cancellation. We computed the Casimir energy for gravity and various 
massless fields, obeying different boundary conditions. In the massive case, we pointed out 
a mechanism for stabilizing compact dimensions, requiring that the massive contribution be 
repulsive, while the massless be attractive. 

We discussed in detail the influence of kinetic terms, localized at the orbifold fixed points, on 
the Casimir energy. We pointed out that these brane-localized kinetic terms can also generate 
stable minima for the radion. In the case of localized kinetic terms of size consistent with NDA 
36fl , one can stabilize the compact dimension at a size several times the fundamental length 



cutoff. For larger brane-localized kinetic terms, that size can be larger. We argued that the 
low-energy perturbative expansion in theories where brane kinetic terms are larger than the 
inverse cutoff scale is consistent (in accord with ref. fl35|). 

Applying our results to the 5-d and 6-d models of electroweak breaking, mentioned above 
E3|, we found that, with the minimal field content, the radion potential generated at one 
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loop due to the Casimir energy is repulsive in both cases. We found that it is possible, by 
adding brane-tension and cosmological constant counterterms to find a stable minimum for the 
radion. 



However, in the supersymmetric 5-d case [22], the sign of the required bulk cosmological 
constant turned out to be positive. Thus, deciding whether this fine-tuning is possible in super- 
gravity would need to wait for a full embedding of the model in 5-d supergravity. Alternatively, 
our result might indicate that a classical mechanism of stabilization is preferred; of course, it 
would also have to be shown to be consistent with 5-d supergravity. 

Finally, we note that, in this paper, we concentrated on static issues, namely the existence 
of stable minima of the radion field(s). Dynamical issues, such as the evolution of the universe 
as a backreaction to the Casimir potentials are of interest as well. It would also be extremely 
interesting to find what radius stabilization mechanism is consistent with supergravity in the 
supersymmetric case. 
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8 Appendix: Evaluation of Casimir sums 

Here, we evaluate some of the sums that appear in the calculation of the Casimir energies in 5- 
and 6-dimensional models, using ^-function regularization. Even though the methods for their 
evaluation are standard, see f|6], |5l|], we show some of the details for completeness. First, we 
evaluate the single sum: 

f(si °' c) '^FWW (78) 

which is convergent for sufficiently large positive Re(s). The idea is to recast the sum into a 
form suitable for analytical continuation to the values of s of interest. Setting a = 0, this is the 
sum that appears in the evaluation of the Casimir energy due to a massive scalar field in 5-d, 
with periodic boundary conditions. For antiperiodic boundary conditions we need a — |. We 
start by noting that eqn. (|7^) defines a periodic function of a with period 1. So we can expand 
eqn. ( |78|) in Fourier series and write: 



F(s;a,c) = V e i2wpa f dye- i27rpy V 
„ Jo „ 
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rn+x 

Y^e^^Y, / dze- i27Tpz 

„ Jn 



n+1 _ . o 1 

[Z 2 + C 2 } s 



/oo I roo _ _ 

r(s) i 

/ rfuu s "fe" u + 2V cos(2vrpa) / 
Jo ~ Jo 



In the third line we used the representation: 



Z - s = — - I dtt s - l e~ zt , (79) 



T(s) Jo 

while in the fourth line we performed the Gaussian integral over z and made the change of 
variable u = c 2 t. In the last line, we recognize the integral representation of the modified 
Bessel function: 

roo 

KJ\x\) = 2 s ~ 1 x~ s duu s ~ l e- {u+ ^ ] (80) 
Jo 

which is valid when Re(x 2 ) > 0. Thus, we finally get for ([781): 



(r (s - i) +4Y(np\c\) s -hos(2npa)K s _ l2 (2np\c\) ] j . 



F(s;a,c) = f^yl c l I 1 ( s ~ ) + 1 L< M<\ ^ c-()s(2^«/)^._ , (27r/;|r| ) | . (81 1 

Note that the expression ([|l]) is now valid even for negative s. According to eqn. (|15|), we need 
to differentiate eqn. ([H]) with respect to s and set s = —2. Since T(— 2) = oo, the derivative 
needs to act only on T(s) in eqn. (jSll). The sum over p, when evaluated at s — —2, can be done 
exactly in terms of the polylogarithm functions 

oo 

U n(x) = ]T — . 

fc=l K 



The result is: 



>7r c| 5 + — (47r 2 |c| 2 Li 3 (g) + 6vr|c|Li 4 (g) + 3Li 5 (g) + h.c.) , (82) 



15 1 1 2vr 4 

s=— 2 



where we used T = — and defined g = e 27 "( a + l l c l). Setting a = or a = ~ (for periodic 
and antiperiodic boundary conditions respectively) and replacing in eqn. fll5|) gives eqn. ( |19"| ) 
of Section ^| (the first term in eqn. fl32"|) corresponds to a finite contribution to the bulk 
cosmological constant and was omitted in eqn. ( |2.4j ), since other infinite contributions have 
already been discarded by the regularization procedure). 
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Next, we turn to the double sum needed in the 6-d models, discussed in Section |j: 

E'l (83) 

where the prime indicates that the zero mode is to be excluded and r = T\ + rr 2 . Writing the 
sum as: 

5 ' ? \n + mr\^ + £ V = E ' E [(„ + mTl )2 + m 2^ + 2 C ( 2a ) ' ( 84 ) 

we note that the sum over n in the first term is of the form eqn. ([78]), with a = rar\ and 
c = mT2- Thus, using eqn. (pl|) we obtain: 

E / , / , 2 . = 2C(2,) + V ) 2) \r 2 \^2t(2s-l) + 
'n + mrr s 1 (s) 



m,n 



+ fTy, N 5 ~* E E {^) co S (27rpmr 1 )^_i(27rpm|r 2 |) . (85) 



Differentiating with respect to s, setting s = —2 and expressing the sum over p in terms of the 
polylogarithm functions as we did in eqn. (|8"2"|), gives eqn. QB2"D of Section [|. 
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